Scalar functions for wave extraction in numerical relativity 
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Wave extraction plays a fundamental role in the binary black hole simulations currently performed 
in numerical relativity. Having a well defined procedure for wave extraction, which matches sim- 
plicity with efficiency, is critical especially when comparing waveforms from different simulations. 
Recently, progress has been made in defining a general technique which uses Weyl scalars to extract 
the gravitational wave signal, through the introduction of the quasi- Kinnersley tetrad. This proce- 
dure has been used successfully in current numerical simulations; however, it involves complicated 
calculations. The work in this paper simplifies the procedure by showing that the choice of the 
quasi- Kinnersley tetrad is reduced to the choice of the time-like vector used to create it. The space- 
like vectors needed to complete the tetrad are then easily identified, and it is possible to write the 
expression for the Weyl scalars in the right tetrad, as simple functions of the electric and magnetic 
parts of the Weyl tensor. 

PACS numbers: 04.25.Dm, 04.30.Db, 04.70.Bw, 95.30.Sf, 97.60.Lf 



I. INTRODUCTION 

Weyl scalars are promising tools for wave extraction in 
numerical relativity, and are being used for this purpose 
in current binary black hole simulations. As has been 
shown in 0, 0, 0, Q , the critical issue is choosing the 
right null tetrad in which to calculate the Weyl scalars. 
The notion of a quasi- Kinnersley frame was introduced, 
which leads to a robust technique for their calculation. 
However, in spite of the relevant result, the procedure 
is complicated in practice, although it has been imple- 
mented with success in @, @, H, [l(| HH, EH • This paper 
gives a fully well-defined and simplified procedure, whose 
results are completely analytical and require no further 
numerical calculations. The pivotal step is the identifi- 
cation of a "preferred" time-like observer. Then, one can 
find simple and well-defined expressions for the scalars in 
the quasi- Kinnersley frame. 



II. WEYL SCALARS 



Weyl scalars are defined as 



*o = C abcd tm b i c m d , 

*i = C abcd e a n b £ c m d , 

* 2 = C abcd £ a m b m c n d , 

*3 = C abc d a n h m c n d , 

*4 = C abcd n a fh b n c fh d , 



(la) 
(lb) 
(1c) 
(Id) 
(le) 



where C a bcd is the Weyl tensor, and £ a , n a , m a and rh a 
comprise a null tetrad with two real and two complex 
conjugate null vectors, such that £ a n a = — 1 and m a rh a = 
1. 

The correct Weyl scalar calculation for wave extrac- 
tion is the one performed in the null tetrad which con- 
verges to the Kinnersley tetrad [HI, [l4[ when the space- 



time approaches a Petrov type D space-time (a particu- 
lar algebraic class of space-times to which Schwarzschild 
and Kerr belong). This tetrad has been dubbed a quasi- 
Kinnersley tetrad. Refs 0, 0, H, 0, @ show that the quasi- 
Kinnersley tetrad belongs to a group of tetrads called the 
quasi Kinnersley frame, whose tetrads are connected to 
each other by spin/boost (type III) transformations. One 
possible quasi- Kinnersley frame was found to be one of 
the three transverse frames where V^i = ^3 = 0. 

The existence of transverse frames has already been 
studied in [l[ and in particular it has been shown that for 
an algebraically general space-time (Petrov type I) three 
transverse frames always exist, while for algebraically 
special space-times (Petrov types II,D,III,N) the number 
of transverse frames varies depending on the algebraic 
type. Hereafter we will however assume to be dealing 
with Petrov type I space-times, which are the ones nor- 
mally considered in a numerical simulation. 

In this article, Weyl scalars are expressed within the 
context of the standard 3+1 decomposition of Einstein's 
equations. A space-like foliation of the space-time is in- 
troduced, with N a the time-like unit normal to the space- 
like hypersurfaces. On each hypersurface, we introduce 
a set of three orthonormal vectors, referred to as u a , x a 
and y a . The set of four vectors - one time-like, the other 
three space-like - leads to a straightforward definition of 
the null tetrad: 



r 

n a 
m a 



2-5 (N a - u a ) . 
2"! (N a + u a ). 
2-i(x a -iy a ), 
2"2 (x a +iy a ). 



(2a) 
(2b) 
(2c) 
(2d) 



The Cauchy foliation of the space-time allows us to 
define the electric and magnetic parts of the Weyl tensor 



2 



E ac — -CabcdN N , 

Bac = —-^(-ab mn C mnc dN b N d , 



(3a) 
(3b) 



where e a b mn is the four dimensional Levi-Civita tensor. 
The final expression for the Weyl scalars is then given by 



vj/ 2 

*4 



- (E bc - iB bc ) m m°, 
2~* {E bc -iB bc )m b u c , 
-2- 1 (E bc -iB bc )u b u c , 
-2"5 (E bc -iB bc )m b u c 

- (E bc - iB bc )fh b fh c . 



(4a) 
(4b) 
(4c) 
(4d) 
(4e) 



Eq. ([4]) is our starting point for determining the correct 
null tetrad for computing the Weyl scalars. Specifically, 
our quasi- Kinnersley tetrad procedure will identify the 
optimal u a , x a and y a for wave extraction, assuming the 
right N a time-like vector has been chosen. 



to calculate them. On the other hand, we know that for 
a specific choice of N a , this condition must hold, as we 
know that transverse frames where = $3 = do in- 
deed exist. This means that there is a "preferred" time- 
like observer who sees those eigenvectors aligned. We do 
not know yet how to identify unequivocally such an ob- 
server, but it is clear that this is the observer to choose 
in order to have meaningful physical results. 

Assuming to have chosen the right time-like vec- 
tor for the calculation of E ab and B ab , and denoting 
{E u ,E x ,E y } the eigenvalues of the electric part, and 
{B u , B x , By} the eigenvalues of the magnetic part, we 
easily obtain (using Eq. (ple|0 the expression for ^4: 



E^ 



. Bit — £> v 



A)TF 



(8) 



Eq. ||5J) indicates that to compute ^4 in the transverse 
frame, one only needs to know the eigenvalues of the 
electric and magnetic parts of the Weyl tensor. 



III. THE TRANSVERSE FRAMES 



IV. EIGENVALUES OF THE ELECTRIC AND 
MAGNETIC PARTS 



First, we calculate the frame where = ^3 = 0. 
Using Eq. (|4b|) and (|4d[) it can be easily shown that this 
condition corresponds to 



E bc u h x c = 0, E bc u b y c = 0, 
B bc u b x c = 0, B bc u b y c = 0. 



(5a) 
(5b) 



The four equations written in Eq. ([5]) can be reduced to 
two by noticing that a spin transformation rotates the 
vectors x a and y a in their plane, leaving the condition 
^1 =^3 = unaltered. Thus, an equivalent way of 
writing Eq. ([5]) is 



E bc u b w c = 0, B bc u b w c = 0, 



(6) 



In this section, we calculate the expressions for the 
eigenvalues of the electric part of the Weyl tensor. The 
procedure for the magnetic part is identical. The equa- 
tion for the eigenvalues is given by 



El 
3 



where Ej and Ej are defined as 

77> 77 TT<lb 77 77' b 771 C TTi a 

The three eigenvalues can be written as 



(9) 



(10) 



where w a is a generic vector lying in the plane perpendic- 
ular to u a . The only way Eq. (|6|) can hold for any vector 
w a perpendicular to u a is by having both products Eb c u b 
and Bb c u b parallel to u a , such that 



E u = -(£ + £*), 

E x = - (e^£ + e^r 



E v 



e — £ + e —£* 



(11a) 
(lib) 

(11c) 



E br U = E„U r , B h rU b 



B u u c , 



(7) where £ is a complex number to be determined. This 



This means that the vector u a is an eigenvector common 
to the electric and magnetic parts of the Weyl tensor, 
E u and B u being the corresponding eigenvalues. The 
three orthonormal common eigenvectors correspond to 
the three possible transverse frames. 

It is evident that the electric and magnetic compo- 
nents of the Weyl tensor cannot have their eigenvectors 
in common for any choice of the time-like vector 7V a used 



construction guarantees that the sum E u -\- E x -\- Ey van- 
ishes, as is expected from the trace-free property of the 
electric part of the Weyl tensor. It furthermore guaran- 
tees that all three eigenvalues are real as expected from 
the symmetric property of E a b- Using this technique, we 
only need to calculate the real and imaginary parts of £ , 
which constitute the two degrees of freedom of the prob- 
lem. Substituting Eq. (jTTJ) into Eq. @, we obtain the 
following relations: 



3 



Ei = 6££*, 



Ej 



3 £ 3 + (£*f 



(12a) 
(12b) 



Writing the expression for £ in terms of modulus and 
phase, 8 = \£ \ e e , the following expressions are derived 
from Eq. $T2§: 



\£\ = 



- arccos 

- > 



-V6 



E ab E b c E ca 
, [E ab E ab Y 



(13a) 
(13b) 



The equations for B are similar: 



161 = 



6 f 



B ab B ab 



() 



- arccos 
3 



V6 



r> r>b r>ca 



[BabB a 



b]2 



(14a) 
(14b) 



Eq. ©, (dU), H31) and (HD give the expression for the 
real and imaginary part of ^4 in the transverse frames 



(* 4 ) TF = - V3 \£ I sin ( B £ + ^ ) , (15a) 



(*4) 



TP 



V^ISIsin f 6b + -77- ] 



3 / 



(15b) 



where k is an integer that can assume the values 
{ — 1,0, 1}, corresponding to the three different transverse 
frames. For later convenience, the expression of ^ 2 in 
this same frame is also given: 



2fc?r 



(* 2 )? F = n cos I H 
(*2)? F = 

For this particular tetrad choice 'J'o = ^4 



(16a) 
(16b) 



is taken to be that given by the Kinnersley tetrad, whose 
expression in Boyer-Lindquist coordinates is 



N a 



1 

71 



a 2 2E-A a 

~ ' 2S ' 0, n 



(17) 



where tt = £ = r 2 + a 2 cos 2 6< and A 



'2 \ lr. M is the mass and a the rotation parameter of the 
black hole. 

The invariant quantities defined in Eq. (|10p are then 
given by 



E, 



Bi = 



Q M 2 r 2 ( r 2_ 3a 2 cos 2 9 y 

(r 2 + a 2 cos 2 9) 6 
6M 2 a 2 cos 2 9 (3r 2 - a 2 cos 2 
(r 2 + a 2 cos 2 9) 6 



(18a) 



-, (18b) 



while E 1 



E 2 £? 2 - 

and Bj — Substituting Eq. |T 



V6 



into Eq. (fT5|) and (fT4")) gives the following result: 



M (r 3 - 3ra 2 cos 2 9) 

(r 2 + a 2 cos 2 6*) 3 
M (3r 2 a cos 9 — a 3 cos 3 
(r 2 + a 2 cos 2 6*) 3 



(19a) 
(19b) 



and Of; = 8b = 0. This leads to the conclusion, us- 
ing Eq. (fT5|) . that ^4 = in the frame corresponding 
to k = 0. Furthermore, we find the correct value for 
"to = -, — J^- — ttt, which is exactly what we expect in 

£ (r+ia cos 8)- 1 ' J ^ 

the Kinnersley frame. We conclude that the frame cor- 
responding to k = converges to the Kinnersley frame 
in the limit of a single Kerr hole; in other words, it is a 
quasi- Kinnersley frame. 

We emphasize again that in our calculation only the ex- 
pression for the time-like vector N a was given "a priori" ; 
no further assumptions on the other space-like vectors 
needed to create the null tetrad were made. In fact, the 
correct identification of those space-like vectors is implic- 
itly achieved using the expressions for the Weyl scalars 
given in Eq. (fT5)) and HT 



V. THE SINGLE KERR BLACK HOLE LIMIT 



VI. THE WEAK FIELD LIMIT 



Eq. (|15p is the expression for ^ 4 in the three transverse 
frames; the validity of this equation is verified below by 
calculating the expression for a single Kerr black hole 
and obtaining the expected results for the Weyl scalars. 
Boyer-Lindquist coordinates are used for the calculation. 

The expressions for the electric and magnetic compo- 
nents of the Weyl tensor are dependent on the time-like 
N a chosen for their calculation. For now, the correct 7V M 



It is assumed here that the metric is of the form g a b — 
Vab + hah where r\ a b is the flat Minkowski metric and h a b 
is the perturbation. 

In the limit of flat space-time, the N a vector of the 
Kinnersley tetrad given in Eq. (|17[) assumes the value 



N a 



. In relation to this vector, E„h and 



2V2' 2V2' u ' u 

B a b are then given, using the transverse-traceless gauge 
for the metric perturbation, by 



4 



ld 2 h 



2uTT 



Baa — — — 



2 dt 2 
1 ° 

2~dt 2 



d 2 h T J 



Ea 



B §4> ~ 



at 5 



2 



dt 2 



(20a) 
(20b) 



where the hatted indices indicate contraction with the 
angular tetrad vectors. 

The expressions for ^4 and ^2 given in Eq. (fT"5f and 
(|16|) simplify considerably in a perturbed flat space- 
time. Specifically, the contribution to the curvature from 
the background vanishes, which means that 'J 2 must 
be zero. Using Eq. (Ti~6)) , one finds that in this case, 
©e = @B = § (assuming to have chosen k = corre- 
sponding to the transverse frame which is also the quasi- 
Kinnersley frame, as shown in the previous section). Sub- 
stituting this result into the expression for ^4 given in 
Eq. (fl5|) . one concludes that in this particular limit, 



(*4) 



QKF 



E nh E ab 



+ i 



B nh B ab 



(21) 



It is interesting to notice that in the limit of per- 
turbed flat space-time, the expression found for \& 4 in the 
quasi- Kinnersley frame (Eq. (|f 5|l ) naturally converges to 
a Poynting vector-like expression. More generally, the an- 
gles 0£ and 0g seem to play a fundamental role in divid- 
ing the curvature expressed by the quantities E a tE ab and 
B a bB ab between background and gravitational waves. In 
this article, we have in fact studied the two limits where 
these angles are either zero, so that the contribution 
to the curvature is given by the background, or -|, in 
which case only the gravitational radiation contributes 
to the curvature. In a general situation, one expects 
< 9 £ < ~ and < Q B < |, so that both the back- 
ground and the gravitational radiation contribute to the 
curvature. 

The expressions for E a bE ab and B a bB ab are given (ne- 
glecting quadratic and higher order terms) by E a bE ab w 

d 2 hl 



and B ab B ab 



, TT 

dt 2 



leading to the 



the two gravitational wave degrees of freedom expressed 
in the TT gauge. The results found in the previous two 
sections justify the use of Eq. (fT5|) as a valid formula for 
wave extraction: it vanishes for a single black hole space- 
time while it is related to the gravitational wave degrees 
of freedom when dealing with a perturbed space-time. 
VII. THE N a TIME-LIKE VECTOR 
DETERMINATION 



In the previous sections, it has been shown that the 
choice of the correct tetrad, using the quasi- Kinnersley 
frame properties, can be reduced to the correct choice 
of time-like vector used for the calculation of the elec- 
tric and magnetic parts of the Weyl tensor. Once this 
is done, the other vectors needed for the construction of 
the null tetrad are easily identifiable and the expressions 
for the Weyl scalars found using these vectors were given 
(Eq. (TJ11) and (USD). 

The "preferred" time-like vector has been found to see 
the electric and magnetic parts of the Weyl tensor with 
their eigenvectors aligned. Work is still in progress to de- 
termine this vector in a general and robust way; however, 
a preliminary valid choice is that of a time-like vector 
which converges to the right one asymptotically, which 
guarantees that the expressions for the scalars are invari- 
ant at first order in perturbation theory. 

In a numerical simulation, the straightforward choice 
for the time- like vector is that of the normal to the space- 
like Cauchy hypersurface. In the limit of flat space-time, 
such a vector converges to N a = [1, 0, 0, 0] (assuming the 
lapse a — ► 1 and the shift f3 a — > in this limit) which 
happens to be boosted in the flat limit with respect to the 
N a found from the Kinnersley tetrad. The result of this 
is that the factor | in Eq. (f2"2")) disappears. This choice 
of the N a time-like vector can be a good approximation. 
However, a way to correctly choose N a in general must 
be found to rigorously complete this work. 



final expression for ^4: 
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(* 4 ) 



QKF 



1 



dt 2 



dt 2 



(22) 



Eq. (f2"2"|) shows that in the limit of a slightly perturbed 
flat space-time, the expression of 'J 4 is directly related to 
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